Skyrmions in a Doped Antiferromagnet 
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Magnetization and magnetoresistance have been measured in insulating antiferromagnetic 
La2Cuo.g7Lio.03O4 over a wide range of temperatures, magnetic fields, and field orientations. The 
magnetoresistance step associated with a weak ferromagnetic transition exhibits a striking non- 
monotonic temperature dependence, consistent with the presence of skyrmions. 

PACS numbers: 75.47.Lx,72.20.My,75.50.Ee 



A remarkable manifestation of complexity in mag- 
netic systems is the emergence of topologically non- 
trivial arrangements of spins, such as skyrmions 
These are "knots" in an otherwise ordered spin tex- 
ture, which behave as excitations with particlelike prop- 
erties. Skyrmions are stabilized by a magnetic field 
B in several ferromagnetic (FM) metals [2|-|5|, where 
they manifest themselves in the electronic transport 0- 
0, H , or they form a detectable periodic skyrmion lat- 
tice [5|. Skyrmions have been predicted to emerge also 
in the ground state of doped antiferromagnetic (AFM) 
insulators but the identification of such isolated 

skyrmions is a challenge. Neutron scattering, for exam- 
ple, would not be a definitive probe, since skyrmions here 
do not form a lattice, while their possible signatures on 
transport may be screened by the insulating character of 
the carriers. 

In a FM system, a charge carrier with the spin aligned 
to the magnetic background preserves its metallic charac- 
ter with a mass renormalization due to scattering by low- 
energy spin waves. Thus, the more complex spin exci- 
tations associated with topologically nontrivial magnetic 
textures have characteristic signatures in transport. This 
is indeed the case for skyrmions in the polarized quan- 
tum Hall state of a two-dimensional electron gas [3, 111 , 
in colossal magnetoresistance (MR) manganites [3], and 
in the three-dimensional FM MnSi [1, In a doped 
AFM insulator, the description of transport is complex 
already in the topologically trivial AFM ground state, 
since the carriers cannot move without inducing spin-flip 
scattering [1, [lH. If, however, the external B field causes 
a change in transport that depends on the configuration 
of the spin background, the MR becomes a key probe to 
identify signatures of anomalous spin textures. This sit- 
uation can be realized in the AFM insulator La2Cu04, 
lightly hole doped with Li. Here we show that the striking 
nonmonotonic temperature (T) dependence of the MR 
step associated with a i3-induced magnetic transition is 
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FIG. 1: (color online) (a) Orthorhombic crystal structure of 
La2Cu04; Li replaces Cu in CUO2 (ab) planes. The spins on 
the Cu atoms have the AFM order with the large staggered 
component no along b and small FM components L = D x no 
along c (D, DM vector, aligned along a). The weak FM mo- 
ments have staggered order along c at B = 0. (b) This AFM 
background becomes deformed near Li sites, via skyrmion for- 
mation. For simplicity, we show the skyrmion configuration 
of the FM moments only. The circle represents the skyrmion 
core size A. (c) Weak FM moments around the Li site across 
the weak FM transition: due to skyrmion formation, spins 
have a FM order inside a distance of the order of A at B < Be 
and AFM outside it, while the reverse is true for B > Be- 



consistent with the skyrmion formation. 

La2Cu04, the parent material of La2Cui_a;Li2;04, is a 
Mott insulator with a residual AFM coupling between the 
S = 1/2 spins located at the Cu^"'" ions, forming a nearly 



square lattice with a small orthorhombic distortion 12 1 
[Fig-da)]. A small antisymmetric Dzyaloshinskii-Moriya 
(DM) exchange between neighboring moments causes a 
uniform canting of the spins, leading to a weak ferro- 
magnetism per Cu02 plane, along with AFM ordering 
[Fig. [TJa)]. When La2Cu04 is doped with charge carri- 
ers, e.g. via Sr or Li substitution, the carriers frustrate 
the magnetism, leading in both cases to a strong reduc- 
tion of the Neel temperature T/v 13|, llj| . Different types 
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of magnetic textures may be stabilized depending on the 
geometry and character of the dopants. For example, the 
diagonal incommensurate magnetic signatures reported 
by neutron-scattering experiments at very low (insulat- 
ing) doping (x) when Sr^+ replaces La"^"*" on top of 
a 4-Cu plaquette can be explained well by the forma- 
tion of local spin spirals [1, At higher x, the incom- 
mensurate magnetic signatures, now rotated by 45°, have 
been consistently attributed to the emergence of spin and 



charge stripes |17| . The substitution of Li"*" for Cu^"*" in- 
plane, on the other hand, may stabilize the formation 
of skyrmions that are obtained, roughly speaking, 
by reversing the average spins in a finite region of space 
around the dopant [Fig. [Ijb)]. At low x, skyrmions are 
expected to affect only the tails of NMR or neutron line 
shapes [lo| , making any unambiguous interpretation dif- 
ficult, but even this has not been reported so far. Here 
instead we take advantage of the weak FM moments to 
tune the orientation of the AFM background using a uni- 
form B. This gives rise to magnetotransport that is sen- 
sitive to the spin configuration around the dopants, and 
consequently to its nontrivial textures. The resulting MR 
data on AFM La2Cuo.97Lio.03O4 are consistent with the 
presence of local skyrmions. 

A high-quality La2Cui_xLi2,04 (Li-LCO) single crys- 
tal with a nominal x = 0.03 was grown by the traveling- 
solvent floating-zone technique [18|. Magnetization M 
was measured in a standard superconducting quantum 
interference device magnetometer in 0.1 < B{T)< 7 ap- 
plied either parallel or perpendicular to the Cu02 (ab) 
planes. The in-plane resistance R was measured on a 
bar-shaped sample with dimensions 2.2 x 0.57 x 0.41 mm'^ 
using a standard four-terminal ac method (~ 7 Hz) in the 
Ohmic regime. The contacts were made by evaporating 
Au and annealing at 700 °C in air. MR was measured 
by sweeping either B \\ ab or B J- ab a,t constant T in 
the 5-190 K range. The sweep rates were low enough to 
avoid the heating of the sample due to eddy currents. 

In La-based AFM systems, the easy axis for the spins 
is the longest (&) of the two in-plane orthorhombic di- 
rections [19|. The DM vector D is oriented along a, the 
AFM order parameter no along b, so that the weak FM 
moments L = D x no are parallel to c [Fig. IHa)]. One 
signature of the presence of the weak FM moments is a 
peak in the magnetic susceptibility at T^r, such that it 
is more pronounced for B || c than for B || b [2^, 21 1. 
This is indeed observed also in Li-LCO [Fig. [2][a)]. Since 
the sample is twinned, the B direction within the plane 
is not specified, but only twins having i? || & contribute 
to anomalies in M (and resistivity p). is suppressed 
by both B \\ c [Fig. [Sib)] and B \\ ab (not shown). More- 
over, the low-T upturn of M is similar to that in un- 
doped and Sr-doped La2Cu04 (LSCO) and thus it 
is irrelevant to the skyrmion formation. A second signa- 
ture of the weak FM moments is found in M{B) at fixed 
T. Because of the weak interplane AFM coupling, the 
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FIG. 2: (color online) (a) M vs. T for B \\ c and B \\ ab. The 
B \\ ab curve is shifted up by 3.5x10"^ emu/g. (b) A4 vs. T 
for several B || c. The curves at 1 T and 0.1 T are shifted 
down by 1.5x10""^ emu/g and lxlO~^ emu/g, respectively. 
Arrows denote the Neel temperature; Tn{B = 0) ~ 180 K. (c) 
M vs. -B II c at different T. The arrows denote the direction 
of B-sweeps. Insets: Spin configuration for (ci) B = and 
(02) B > Be, where Be is the critical field for the spin-fiop 
transition (red and blue arrows: full Cu spins; green arrows: 
weak FM moments), (d) M vs. B || ab at T = 5 K. Inset: 
Continuous rotation of spins in the be plane for -B || &, where 
Be is the saturation field. Dashed lines guide the eye. 



weak FM moments have a staggered order along the c 
axis [Fig. [2lc) inset (ci)]. A sufficiently large field Be \\ c 
can overcome the interplane AFM coupling and induce 
a discontinuous spin-fiop reorientation in both the weak 
FM and the large AFM components [Fig.[2Ic) inset (C2)], 
causing the so-called weak FM (WFM), first-order transi- 
tion M,IM[2l. This results in a jump AM{T) at B^T) 
22, 2| [FigJ2i;c)l Both AM{T) and BJT) decrease 
with T [20|, I22I, |23|, following the decrease of the stag- 
gered magnetization no(T) due to thermal fluctuations 
For B II b, the weak FM moments induce a con- 
tinuous rotation of no in the be plane [l9|, Hil (Fig. [^d) 
inset), so that M{B) increases smoothly [Fig.[2ld)]. Thus 
the magnetic properties of Li-LCO resemble very much 
those of Sr- and O-doped La2Cu04 [l^, [El, [ij where 
skyrmions cannot form. 

At very low x, the motion of the holes may be described 
as the hopping of (spinless) charges in the two AFM sub- 
lattices [sl, [ll|. Because of orthorhombicity, the majority 
of holes have momenta close to (7r/2, — 7r/2) [i^, i.e. the 
most relevant pocket is the one along the b direction. 
When the Coulomb potential provided by the dopant is 
taken into account, the holes get localized, with a typical 
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FIG. 3: (color online) (a) p vs. T ak. B = 0. (b) MR for 
B II ah aX several T. (c) MR vs. T at 18 T for B || ah. 
Inset: The sketch of the weak FM moments (blue arrows) 
and the corresponding inter layer hopping t± for B \\ ah {to, 

(2) 

maximum value of fx; Be , saturation field), (d) MR for 
B II c at three selected T. Arrows denote the directions of 
the _B-sweeps. Inset: Weak FM moments (blue arrows) before 
and after the weak FM transition at Be and the corresponding 
field-induced tx. (e) MR for B || c at 15 < r(K)< 190. (f) 
The nonmonotonic T dependence of the MR step size shown 
in (e), and the monotonic decrease of the magnetization step 
with T. The dashed line guides the eye. The calculations 
of the MR step size in the presence and in the absence of 
skyrmions §31 are also shown. 



2D pancake-like exponential envelope ^(r) ^ exp(— r/^o) 
(r, planar distance from the impurity site; localiza- 
tion length) The low-T electronic transport is then 
expected to occur via variable-range hopping (VRH) be- 
tween localized states [2^ at a characteristic hopping 
distance rh{T). Indeed, the insulating behavior of R 
[Fig. ma)] follows a 2D VRH form R cx exp(To/T)i/3 
(To — 3753 K) for T < 23 K, and it crosses over to 
R (X e^]i{EA/T) {Ea = 313 K) at higher T similar 
to the behavior of other La-based AFM samples [l^ ■ 

The in-plane MR [R{B)/R{0) - 1] at first glance also 
resembles that of AFM LSCO [H, For example, 

for B II ab, the MR is negative, decreases continuously 



with a tendency to saturate above some threshold field 
[Fig-EIb)], and its overall magnitude increases monoton- 
ically as T is reduced [Fig. EI^c)]. For B \\ c, MR is also 
negative, with a steplike decrease [Fig. El^d)] at the same 
critical field Be where the WFM transition occurs and 
the uniform M shows a jump [Fig.[2Jc)]. Moreover, both 
M [Fig. Wic)] and MR [Fig. [S^d)] exhibit a hysteresis as- 
sociated with a first-order phase transition. However, as 
shown below, it is the nonmonotonic T dependence of the 
MR magnitude [Fig.lHIe)] and \AR\/R [Fig.EJf)], the MR 
step size at Be, that signals the formation of skyrmions 
in the ground state of AFM Li-LCO. At high T, \AR\/R 
grows with decreasing T as expected, but unlike AFM 
LSCO [22I, it exhibits a dramatic reversal of the behav- 
ior below r w T5 = 35 K, and starts to decrease as T is 
reduced further. We attribute this decrease to the pro- 
gressive emergence of skyrmions of increasing core size 
that suppress the interlayer tunneling process responsi- 
ble for the negative MR in doped AFM La2Cu04 [ist - 

In general, holes can only tunnel between planes along 
the direction of FM order of the spins: 6 at i? = 
[Fig. [TJa)] and a above the WFM transition. How- 
ever, since the hole momentum is close to (7r/2,— 7r/2) 
due to orthorhombicity, this hopping process is sup- 
pressed at i? = 0. It becomes possible only above the 
WFM transition [27] , modifying the in-plane wave func- 
tion as ^B{f) = /(j't B)4'oif) ("00,-8 are the wave func- 
tions at zero and finite B, respectively). In the hopping- 
conductivity scheme described above, R ^\ ipirt) |^^, so 
that the MR is given by = _(i_g£i^). At 

the typical length scale rh{T) for hopping, f{rh,B) ~ 1-f 
a {rh t^{B)Y > 1 [i^! leading to a negative MR, as ob- 
served in Figs.|3l^b) and|3l^d)-(e). Here t± is the interlayer 
hopping and a = 2/ (CoCo)^ (eo, the localization energy of 
the hole bound state). When B \\ b, t^{B) = tosin6'(B) 
increases continuously with B (Fig.lSjc) inset); here 9{B) 
is the angle with the b direction that spins form in their 
smooth rotation in the be plane (see also Ref. [lij). How- 
ever, for B II c, t±{B) jumps discontinuously (Fig. ^d) 
inset) from zero to the maximum value at the critical 
field Be. Hence, the MR is continuous |2J] for -B || 6 
[Fig. El^b)], while it is discontinuous 22, l23| for B || c 
[Fig. Eld)]. 

In this picture, the T dependence of the step size 
\AR\/R is determined by rh{T) and t±{T). As the trans- 
port crosses over from the activated to the VRH regime 
(where r,,(r) - Cq(Tq/T)^^^) at low T, rh{T) increases 
with decreasing T, and one expects only a monotonic in- 
crease of the MR magnitude. This is indeed observed 
for B II b [Fig. [Sic)] and in AFM LSCO [H]. Since the 



MR behavior is intimately related to the AFM order, 
any mechanism that leads to a global reduction of the 
AFM order parameter no(T) may be expected to con- 
tribute also to a suppression of \AR\/R. However, the 
low-T drop of \AR\/R cannot be ascribed to any such 
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mechanism, since the magnetization step AM(T), which 
is proportional to no(r) [l^, increases continuously in 
the same T range [Fig.[31[f)]. This rules out, for example, 
spin glassiness, which is known to reduce AM(T) in AFM 
LSCO [23], but which here sets in only below 7-8 K [H. 
It also rules out a structural change of the kind observed 
in Lai.79Euo.2Sro.oiCu04 [1^, although none have been 
reported on Li-LCO in our T range. Charge glassiness, 
which gives rise to a positive MR, is also not relevant, as 
it sets in at much lower T [s^l- Thus the only mecha- 
nism that can lead to a decrease in \/S.R\/R is the local 
suppression of tj_, as it happens when a skyrmion forms. 

Close to the Li impurity position r = 0, where the 
skyrmion is centered, no is nearly fully reversed with re- 
spect to the equilibrium direction and recovers only at 
distances r > A [A, the characteristic length scale fixing 
the core size of the skyrmion, Fig.[IJc)]. At i? = 0, the or- 
dering of spins in the two neighboring layers, the one with 
the impurity and the one above it, is nearly FM within a 
distance of the order of A. For B > a full reversal of 
both weak FM and staggered moments occurs: the order- 
ing of the magnetic moments within a distance A from the 
impurity is AFM [Fig. HJc)] , and t_L is suppressed for dis- 
tances up to r ~ A. Thus, if the formation of skyrmions 
suppresses t±_ in f{rh,B), the MR step size is expected to 
decrease. By modeling the gradual increase of A as T is 
lowered [27|, we have reproduced [Fig.[3{f)] the observed 
nonmonotonic behavior in \ AR\/R for B \\ c, character- 
ized by a pronounced downturn below Ts — 35 K. This 
behavior is contrasted to the smooth and monotonic in- 
crease in the case where no skyrmions are formed [25| . 
Since A obtained in our analysis is of the order of one 
lattice spacing [531, such nearly pointlike skyrmions will 
interact very weakly among themselves and they are not 
expected to cause significant changes in bulk quantities, 
such as M . They will be crucial, however, in suppressing 
local processes, such as the vertical tunneling around the 
impurity, associated with the negative MR effect. 

Skyrmions carry a nonzero topological charge, Q ^ 0, 
and a skyrmion configuration is orthogonal to the Q ^ 
Neel state. In order to stabilize a skyrmion at T = 0, the 
hole state must modify the topology of the Cu02 lay- 
ers. Since the spins on the four nearest neighbors to the 
Li sites belong to the same AFM sublattice [Fig. [DJa)], 
the localized hole wave function can be labeled by the 
eigenvalues of the orbital angular momentum: £ = ±1 
(excited states) or ^ = ±i (lowest energy states) [§| . The 
£ — zLi eigenstates have circulation associated with them, 
describing a localized hole orbiting around the Li center 
0. The associated, nonzero spin current induces the 
skyrmion deformation of the AFM background [^, [lo| . 
This circulation might, or might not, be affected by an 
applied B, depending on whether the original symmetry 
of the problem is preserved. For i? || c, the symmetry 
is preserved, £ = ±i are still good quantum numbers, 
and skyrmions survive the applied B. For B \\ ah, on the 



other hand, the symmetry is broken, the £ = ±i states be- 
come mixed, causing the quenching of the hole angular 
momentum, and eventually skyrmion formation is sup- 
pressed. As a result, \AR\/R increases monotonically as 
T is lowered [Fig. [3^c)], as expected without skyrmion 
formation [see theory curve in Fig. El^f)]. Such a mono- 
tonic increase is indeed observed in both B \\ ah and i? || c 
for Sr doping [2^ . where skyrmions cannot occur since 
the hole orbital angular momentum is always quenched, 
from the start, by the mixing between the two sublattices 
around the position of the Sr dopant. For vanishingly 
small B II ah, however, skyrmions might still be favor- 
able and induce a nonmonotonic behavior for |Ai?|/i?, 
but probably for T much lower than those accessed in 
our experiment. Thus, what matters for the skyrmion 
formation is the local structure of the hole wave function 
at the Sr or Li sites, not the way transport occurs. In- 
deed, in the AFM phase of interest, both LSCO [23] and 
Li-LCO show similar VRH or activated behavior, with 



differences emerging only at higher doping 13|, [14 1 . 

In summary, the low-T magnetic and transport prop- 
erties of the AFM La2Cui_j:Lij:04 provide the first ex- 
perimental support for the predictions of skyrmions in 
AFM insulators. Our work may offer new insights into 
the mechanisms that can stabilize or suppress topological 
excitations in complex magnetic systems. 
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The insulating behavior of the in-plane resistance 
R of the antiferromagnetic La2Cuo.97Lio.03O4 sample 
[Fig. 3(a)] follows a 2D VRH form R oc exp(ro/r)i/3 
for T < 23 K (Fig. SI left). At higher T, it crosses over 
to an activated law R cx exp(£'^/T) (Fig. SI right). 
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T-i/3(,(-i/3j t\10"'K-') 

FIG. SI (color online) Left: 2D VRH behavior of p(J) for 
1 < r(K)< 23. Right: Activated behavior of p{T) for 80 < 
r(K)< 190. 

Figure S2 shows the weak FM transition in the mag- 
netization M m B \\ c [Fig. 2(c)] at several temperatures 
and in magnetic fields all the way up to 7 T. 
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FIG. S2 (color online) M vs. i? || c at several T. The arrows 
denote the direction of B-sweeps. 

Here we discuss in more detail the theoretical calcula- 
tion of the magnetoresistance step size as a function of T 



[Fig. 3(f)]. The calculation is based on the known mech- 
anism of the negative MR in doped AF La2Cu04 [1]. By 
using a phenomenological model to account for the pres- 
ence of skyrmions in the case of Li doping, we show that 
the T-dependence of the MR step size is reversed at low 
T, in agreement with the data [Fig. 3(f)]. 

It is widely accepted that the magnetic properties of 
AF undoped La2Cu04, including the role of DM interac- 
tion [2] , can be fully accounted for by a semiclassical con- 
tinuum field theory (the quantum nonlinear sigma model) 
for the local spin moments on the Cu [3]. At very low 
doping, where the static and long-ranged AF order is still 
dominant, a semiclassical treatment of the hole motion 
and of the magnetic textures induced by the carriers is 
still justified. Within this approach, it has been argued 
that the geometry and character of the dopants can in- 
duce different local distortions of the AF background, in 
agreement with the experimental observations. For ex- 
ample, the formation of periodic spin spirals [4], which 
is also supported via cluster Monte Carlo method [5], 
can explain the incommensurate neutron response [6] in 
a Sr-doped material very well. On the other hand, the 
substitution of Li+ for Cu^''" can, according to exact di- 
agonalization studies [7] , stabilize a skyrmion, leading to 
the precise signatures in magnetotransport discussed be- 
low. Finally, at higher doping, when long-range magnetic 
order is completely suppressed and the charge transport 
has a metallic character, various short-ranged spin and 
charge modulations may emerge [8]. 

Here we focus on the case of very low doping. To 
evaluate the MR in the presence of skyrmion defects, 
we extend to the case of Li doping the approach de- 
veloped in Ref. [1], which has been proved to explain 
both qualitatively and quantitatively the MR data in Sr- 
doped AF La2Cu04 [10 ,11]. As a starting point, we 
treat the motion of the hole in the 2D AF background 
within a semiclassical approach to the t — t' — J model [5], 
where the hole spin is constrained to be aligned to the 
AF background, so that the hole motion can be reduced 
to the hopping of a spinless charge in a given sublat- 
tice, with hole momentum near either of the two pock- 
ets (7r/2, ±7r/2). The sublattice degeneracy is further re- 
moved by orthorhombicity of the crystal, which is large. 
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and we assume that only the b pocket {tt/2, —tt/2) is 
occupied [9]. At low doping, the trapping Coulomb po- 
tential V{r) provided by the dopants (either Sr or Li) is 
unscreened, and the transport has an insulating charac- 
ter. The spin degrees of freedom are described using the 
quantum nonlinear sigma model [3] in the presence of the 
DM interaction [2] , which successfully describes the evo- 
lution of the magnetic order as a function of temperature 
T and the applied magnetic field B. 

The 2D picture for the hole motion is partly modi- 
fied by the presence of B that tends to align the weak 
FM moments and, as a consequence, induces free in- 
terlayer tunneling t±{B) of the holes. Indeed, as dis- 
cussed in detail in Ref. [1], hopping between planes is 
constrained by the direction of ferromagnetic alignment 
of the spins. At B = 0, the spins in neighboring lay- 
ers align antiferromagnetically within the ac-plane and 
ferromagnetically within the &c-plane, while the oppo- 
site occurs above the spin-flop transition, induced by the 
alignment of the WF moments along the applied field. 
As a consequence, the hole dispersion is corrected as 
(5ek — —t± cos kz cos{kx ± ky)/2, with the minus sign at 
B < Be and the plus sign at B > Be- However, since 
only the b pocket is occupied, the hole momentum is near 
(7r/2, — 7r/2), leading to a vanishing out-of-plane hopping 
for B < Be and to a finite one for B > Be- The resulting 
Schrodinger equation for the hole wave function can be 
written as a layered model [1]: 

~ v{r)6no) r-t^iB)[r+'+r-'~2yjn] = e^, 

\ 2m* J 

(1) 

Here -0" is the envelope wave-function of the hole in the 
n-th plane, and m* ~ 2me is the mass of low-energy 
quasiparticle excitations in the AF background. For a 
shallow level the exact form of the local potential is 
not important, so we can use a delta-like approxima- 
tion V{r) = —gS{r) for the localizing potential V{r) 
provided by the Li dopant that acts only in the n = 
plane where the hole resides. For t± = 0, the wave func- 
tion of the hole with the full 2D kinetic operator is 
■0o(?") ~ ^o(''/Co)j where Kq is the modified Bessel func- 
tion of the second kind and is the localization length, 
determined by m* and g. For any practical purpose, what 
matters to the magnetotransport mechanism [1] is only 
the asymptotic exponential decay of the hole wave func- 
tion, ^jjo{r) ~ e~'^^^° at r > ^0- Thus, to allow for an ana- 
lytical treatment of the skyrmion case, we shall consider 
in what follows the reduced one-dimensional problem as- 
sociated with the radial part of the envelope function. 
The strength g of the localizing potential is adjusted to 
reproduce for t± — the localization length = 15 A 
extracted from the measured Tq. In particular, Tq was 
obtained from the fit in Fig. 3(a) inset (ai) and was 
estimated from the expression = [13.8/(fcToA^2)]^^^ 
vahd for 2D VRH [12], with the 2D density of states at 



the chemical potential N2 ~ (0.7-4.5) x 10^^(cm^ eV)"^ 
{N = (2-13)/(eV cell) [13]). 

As we explained above, the t± term in Eq. ([T} accounts 
for the interlayer hopping process that becomes possible 
when the in-plane or out-of-plane magnetic field is ap- 
plied, causing a reorientation of the spins in neighboring 
planes. In the perpendicular-field geometry and without 
skyrmions, t±{B) = for B < Be, while t±{B) = to for 
B > Be, where Be is the critical field for the spin-flop 
transition. We note that this description of the hole mo- 
tion is effective only in the presence of long-ranged AF 
correlations, so that MR is expected to vanish at T > T/v. 
While a full treatment of the MR at high T is beyond the 
scope of the present work, we partly accounted for this 
effect by rescaling t± with the AF order parameter no(T), 
which has been calculated following Ref. [2]. Note that 
this does not affect the MR temperature dependence in 
the regime where skyrmions form, because no(T) here 
has a weak T dependence, in agreement with AM{T) cx 
no(r) [2] shown in Fig. 3(f). 

When skyrmions are not present, one recovers the case 
discussed in Ref. [1] for Sr-dopcd La2Cu04. Within a 
VRH scheme, the decrease of the resistance across the 
spin-flop transition is evaluated through an increase in 
the hole's probability |^/'B(r)|^ to propagate at the VRH 
distance = (^o/3)(ro/T)i/3, hence = -(1- 

TT^^^Hw). We note that even outside the regime where 
the resistivity has strictly VRH character {i.e. above 
23 K in our case) what matters for the T-dependence of 
the MR is just the increase of the hopping probability 
r/j(r) with decreasing T. Since we are interested in dis- 
cussing the role of skyrmions below ^ 35 K, we will use 
for simplicity the VRH expression for rh{T) at all tem- 
peratures. This accounts qualitatively for the MR also at 
higher T, where a quantitative agreement with the data 
is not expected within the present approach, as explained 
above. 

The formation of the skyrmion defect around the 
dopant site r — Q affects this interlayer hopping pro- 
cess in a crucial way. Indeed, in the skyrmion configu- 
ration the spins are reversed at the impurity site r = 0. 
Therefore, before the flop, spins in neighboring planes are 
ordered ferromagnetically inside a distance of the order 
of the skyrmion core size A, while the opposite occurs 
above the spin- flop transition [see Fig. 1(c)]. This can 
be modeled by means of a distance-dependent tunnel- 
ing process, i.e. t±{r,B) — to9{r — X) at B < Be, and 
t±{r,B) = toO{X -r) at B > Be, where 9{x) = 1 for 
X < and zero otherwise, with to ~ 3.4 meV [1]. Within 
this simplified, hard skyrmion model for the interlayer 
hopping, according to which the spin quantization axis 
is reversed discontinuously across the core size and the 
skyrmion tail is removed, the parameter A turns out nu- 
merically much smaller than the typical values obtained 
in Ref. [7], where A w 2 — 3 in units of the lattice spac- 
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ing. The above hard skyrmion model could be refined, for 
example, by considering explicitly the smoother r depen- 
dence of a soft skyrmion in the spin-dependent interlayer 
hopping, t±_ = (tg / ^/2S) ^y S'^ + Si ■ S^+i , with the spin 
quantization axis for spins at the «-th layer, S^, being 
reversed continuously across the core size. In this case, 
and as considered in [7] , one naturally expects larger nu- 
merical values for A. It is important to emphasize, how- 
ever, that if one includes also all magnetic anisotropies, 
the antisymmetric Dzyaloshinskii-Moriya and XY inter- 
actions, as well as the interlayer exchange Jj_, the energy 
of the skyrmion far field will be increased and the core 
size will be automatically reduced, as shown in [7]. In 
other words, since the magnetic anisotropies enforce the 
6-orthorhombic direction as the easy axis for the stag- 
gered moments, the core size is automatically reduced in 
order to minimize the energy to be paid, by the circulat- 
ing hole, to generate a skyrmion distortion of the back- 
ground, thus justifying our choice for the hard skyrmion 
description. In any case, since both models lead to anal- 
ogous qualitative physics differing only in the numeri- 
cal value for A, which in both cases is smaller than the 
localization length we shall adopt, in what follows, 
the hard skyrmion model, for the sake of simplicity. We 
also note that a small /lard-skyrmion core size is expected 
from the T dependence of the magnetization step, which 
is not affected by the skyrmion formation. Indeed, a sig- 
nificant deformation of the spin background due to large 
skyrmions would decrease also the value of the AF or- 
der parameter uq , proportional to the magnetization step 
shown in Fig. 3(f). At the same time, a point-like size 
of the magnetic defects justifies also the neglect of the 
interactions between them, as it is implicitly done here. 



By solving again the model [Eq. ([T])], one can see 
that the size of the MR step shrinks due to the pres- 
ence of skyrmions, which effectively reduce the inter- 
layer hopping channel opened by the spin reorientation. 
The size of this effect depends on the size A of the 
skyrmion core, which we expect to increase with de- 
creasing T to lower the energy of the quantum defect 
with respect to its classical value [14]. In Fig. 3(f), we 
modeled the gradual increase of A with decreasing T as 
A(r) = Ao/(l + exp((r - T,fe)/r„), where Ao = O.OS^o, 
Tsk = 17.5 K and = 15 K. Here Tsk represents 
the temperature scale where long-wavelength transverse 
fiuctuations of the AF order parameter (that are con- 
trolled by magnon gaps around 10-20 K) disorder the 
skyrmion ground state. The resulting T-dependence of 
the MR step size is in an excellent agreement with the 
data [Fig. 3(f)]. 
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